theta1 = rand(3,1)
theta2 = [0.1 0 0]';

vee(logm(expm(so3(theta1)) * expm(so3(theta2))))

theta1 + (jacobian_right(theta1))\theta2
